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This report summarizes recent developments in the analysis of various waveguide com- 
ponents and feedhorns using Modal Analysis (Mode Matching Method). A brief descrip- 
tion of the theory is presented , and the important features of the method are pointed out. 
Specific examples in circular , rectangular, and coaxial waveguides are included, with com- 
parisons between the theory and experimental measurements. Extensions to the methods 
are described. 


I. Introduction 

Modal Analysis has been shown to be a highly accurate and 
versatile method for analyzing a wide variety of waveguide 
devices [l]-[4]. The method is capable of accounting for 
multiple reflections within the device, stored energy at each 
discontinuity, and higher-order mode propagation if it occurs. 
Its high accuracy makes it useful for tolerance studies after a 
final design has been determined. This report compares com- 
puted and experimental results for the scattering parameters 
of three examples. One example is taken from each of the 
following waveguide types: rectangular waveguide, circular 
waveguide, and coaxial waveguide propagating the TE n mode. 

II. Theory 

The theory described below is well known and is summa- 
rized in [1] . In applying the modal analysis method, the wave- 
guide device is broken up into a series of sections that are 
joined by a step discontinuity as is shown in Fig. 1 . For smooth 
changes in waveguide dimensions, the change is approximated 


by a large number of steps. At this point, the type of wave- 
guide is arbitrary but the common area between the two 
guides must be identical to the cross-section of the smaller 
waveguide. This eliminates a class of offset connections but 
is usually not important for analyzing a practical device. In 
addition, for the circular waveguide and the coax, all guides 
are required to possess the same center line. This simplifies 
the analysis since only modes with one azimuthal variation 
need to be considered. Again, this is not restrictive for most 
practical applications. Next, the development of some of the 
important equations is presented. In Fig. 1, the fields to the 
left of the junction (z < 0) are represented as a sum of the 
normal modes of waveguide I. 
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Here M is chosen large enough for convergence, and e Im 
and h Im are the normalized vector functions for the mth mode. 
For example, in a circular waveguide, m = 1 = TE n , m = 2 
= TM n , m = 3 = TE 12 , etc. A Im represents the magnitude of 
the forward traveling mth mode, and B Im the magnitude of 
the reverse traveling mth mode. 


i 


The normalization of e Im and h Im is such that 

(?/m X — /m^ * ds = R mm 
and from the orthogonality of the waveguide mod 
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Similarly, in region II 


Un 


N 


-II * ^Iln e 




'Un 


Un 


( 5 ) 


( 6 ) 


where N is the number of modes chosen in region II. Rela- 
tions analogous to Eqs. (3) and (4) hold for this region also. 

Matching the electric and magnetic fields over the common 
aperture results in the following scattering matrix equation. 
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In Eqs. (8) and (9), B Y and B n are vectors containing the 
reflected-mode amplitudes, while ^4j and A u contain the 
incident -mode amplitudes. The derivation of these equations 
is given in the Appendix. 

For the normalized vectors, the power carried by the mth 
forward traveling mode is given by \A Im | 2 , and for the mth 
reverse traveling mode is I B lm | 2 . 

Next, the matrices for a straight section of length L are 
needed. The solution is trivial, giving 


[S n ] - [0] 

(ID 

[S 12 ] - 

(12) 

P ai l = M 

(13) 

[S 22 ] = [0] 

(14) 

where [0] is the zero matrix and [y l2 ] and [ 7 2] ] 
matrices with elements 

are diagonal 

y 2l (n,n) = e y " L = y 12 (n,n) 

(15) 


y n being the propagation constant for the nth mode and L 
being the section length. 

This completes the summary of the required equations for 
each step in the analysis. Using these results, matrices for each 
step and each straight section in the device are determined and 
then combined using equations in [1] . At the completion of 
the analysis, the overall matrix is obtained, relating the nor- 
malized output vectors B l and B n at the ends of the device to 
the normalized input vectors A x and^4 n . 
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In many situations, a set of modes is incident only on the 
(9) left end of a device and one wants to determine the reflected 
and transmitted modes. The user specifies the input mode 
vector ,4j, A u = 0, and Eqs. (16) and (17) become 
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III. Results 

Computer programs [5] have been written to carry out the 
above calculations for rectangular, circular, and coaxial wave- 
guides. For the junctions involved, the integrals of Eq. (A-6) 
can be carried out in closed form, which greatly simplifies the 
programming. Three examples, one in each waveguide type, 
are used to demonstrate the excellent agreement between theory 
and experiment. All of the measurements were made using an 
HP 8510 network analyzer and a full two-port calibration. 

The first example, shown in Fig. 2, is a circular waveguide 
transition. For the theoretical results, 25 modes were used in 
the input section. The number of modes used in subsequent 
sections is chosen by the program for optimum convergence 
[1 ] . The return loss measurements (Fig. 3a) are typically within 
0.2 dB, except near the minimum reflection point at 8.25 GHz. 
Phase results (Fig. 3b) are also in close agreement, typically 
within a few degrees across the band. For nearly every obser- 
vation point, the difference between theory and measurement 
is within the accuracy specification of the network analyzer. 
Slight inaccuracies in the waveguide dimensions and rounding 
of some of the corners can also account for the small disagree- 
ment that remains. Figure 4 illustrates the convergence of the 
solution at 8 GHz as the number of modes used in the input 
section is increased. From these plots, we see that the solution 
has stabilized once 20 or more modes are used in the input 
waveguide. The number of modes required for convergence 
depends on the particular device, but in general larger wave- 
guides with respect to a wavelength and thin irises require that 
more modes be used in order to get the same accuracy. 

A rectangular waveguide example is shown in Fig. 5 and 
both theoretical and experimental results are given in Fig. 6. 
The device is a WR125 to 0.8-inch square-to-WR125 transition 
that was fabricated for use in the ring resonator at 8.51 GHz. 
The device consists of nine waveguide sections. The figure 
shows that, as in the previous example, the theoretical and 
experimental results are in excellent agreement. Only slight 
discrepancies appear near the minimum reflection point. More 
modes may be needed to represent the field in this region, or 
else the 0.030-inch radius on all corners, which was not 


accounted for in the calculation, may have a stronger effect 
in this frequency band. For this example, modes with first 
index m less than or equal to 7 and second index n less than 
or equal to 6 were used in the input guide. As with the circular 
waveguide program, maximum mode indices in the following 
sections are chosen according to waveguide size and symmetry 
considerations. 

The final example is the coaxial iris shown in Fig. 7, with 
theoretical and experimental return loss results shown in 
Fig. 8. The coaxial region is excited by a TE n circular wave- 
guide mode that excites only the higher-order coax modes 
with first index equal to 1 ; the normal TEM coax mode is not 
excited in this case. Measurements of the iris were made by 
calibrating in the circular waveguide and using the time domain 
gating features of the HP 8510 network analyzer to isolate 
the reflections from only the iris. The only other complica- 
tion associated with the coax is that a transcendental equa- 
tion must be solved for each mode in each section in order 
to determine a cutoff wavelength. This increases the computa- 
tion time required to solve a coax problem compared to a 
similar circular or rectangular waveguide problem. As with 
the previous examples, the agreement between theory and 
experiment is good, particularly considering the errors intro- 
duced by using the time domain features of the HP 8510. 

IV. Conclusion 

Three representative examples have been given to demon- 
strate the accuracy of the modal analysis method. A large 
number of waveguide devices such as horns, corrugated wave- 
guides, transitions, filters, and smooth tapers can be analyzed 
using these programs. In addition, several extensions have 
been made to the codes in order to allow for differing dielec- 
tric constants in the sections, making them useful for window 
design. For large smooth-wall or corrugated horns, the reflec- 
tion at the aperture may be neglected, and the far-field pattern 
can be found from the propagating modes in the aperture. In 
addition, the important case of ring-loaded slots [6] , which is 
a combination of the coaxial and circular program, has also 
been programmed, but no experimental results are presently 
available. 
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Fig. 4. Convergence for circular waveguide example: (a) return 
loss and (b) phase. 



Fig. 5. Rectangular waveguide example. 



Fig. 6. Rectangular waveguide return loss results. 
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Appendix 

Derivation of the Waveguide Scattering Matrix Equation 


To derive Eq. (10), the electric fields inside the common 
aperture between the two regions are matched. 

Ej = E u inside Sj (A-l) 

Ej X h Un = E u X h Un inside S 7 (A-2) 



ffi.x 



®„X *„.)•* 


Un* 


(A-3) 


Since E = 0 on the conductor making up the surface S u - S v 
the integral on the right-hand side of Eq. (A-3) may be extended 
over S jj. 
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Equations (A-5) and (A-10) may be recast into a more com- 
pact matrix form, giving 


[P] = [Q] (A n +B n ) (A-l 1) 

[R] 04,-i?,) = [/>] T (£„-/!„) (A- 12) 


Here [P] T is the transpose of the matrix [P] , and [R] is an 
mXm diagonal matrix, and [Q] is an n X n diagonal matrix. 


Using the properties in Eqs. (l)-(6) the following is obtained: 
M 

V (A T +B r )P = G4 tt + P tt )Q (A-5) 

/ , v Im Im' m n v Iln lln'^nn v 
m- 1 

where 


Next, Eq. (A-12) is converted into a scattering matrix for- 
mat relating the normalized output vectors B Y and B n to the 
normalized input vectors A x andvt n . 

The submatrices [S n ] , [S 12 ] , [S 2l ] , and [S 22 ] are derived 
from the [P] , [P] T [P], and [ Q ] matrices by simple matrix 
math and Eqs. (A-l 1) and (A-12). 
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(A-6) is n ] = a*] + fi t Mr 1 m - [/ > i T t^D 

(A-l 3) 

[S i2 ] =2K/R] ([R] + [P] T [R]) [R] T [VS ]" 1 (A- 14) 

(A- 7 ) 

P 21 ] =2 [VS] ([Q] + [P] [P] r ) [P] [VR ]" 1 (A- 15) 


The other boundary condition needed is 
- H u within S 2 


(A-8) 


[S 22 ] = IQ] (IQ] + [P] [P ] 7 )' 1 ([Cl - [R] [R] t ) [VS ]- 1 

(A- 16) 


Following a similar line of reasoning 



(A-9) 


In these equations, [/] represents the unit matrix and 

[VR] [VR] = [R], and [VS] [VS] = [VS] are from 
Eqs. (3) and (7). These factors form the normalization of the 
vectors A and B. This completes the solution for the junction 
between the two different waveguides. 
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